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Observational Strategies of CMB Temperature and Polarization 
Interferometry Experiments 
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ABSTRACT 

We have simulated the interferometric observation of the Cosmic Microwave Back- 
ground (CMB) temperature and polarization fluctuations. We have constructed data 
pipelines from the time-ordered raw visibility samples to the CMB power spectra which 
utilize the methods of data compression, maximum likelihood analysis, and optimal 
subspace filtering. They are customized for three observational strategies, such as the 
single pointing, the mosaicking, and the drift-scanning. For each strategy, derived are 
the optimal strategy parameters that yield band power estimates with minimum uncer- 
tainty. The results are general and can be applied to any close-packed array on a single 
platform such as the CBI and the forthcoming AMiBA experiments. 

We have also studied the effect of rotation of the array platform on the band power 
correlation by simulating the CBI single pointing observation. It is found that the 
band power anti-correlations can be reduced by rotating the platform and thus densely 
sampling the visibility plane. This enables us to increase the resolution of the power 
spectrum in the ^-space down to the limit of the sampling theorem (A£ = 226 ~ tt/^), 
which is narrower by a factor of about \/2 than the resolution limit (A^ ~ 300) used 
in the recent CBI single pointing observation. The validity of this idea is demonstrated 
for a two-element interferometer that samples visibilities uniformly in the uu-annulus. 

From the fact that the visibilities are the Fourier modes of the CMB field convolved 
with the beam, a fast unbiased estimator (FUE) of the CMB power spectra is developed 
and tested. It is shown that the FUE gives results very close to those from the quadratic 
estimator method without requiring large computer resources even though uncertainties 
in the results increase. 
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1. Introduction 

It is theoretically expected that the CMB is polarized. The CMB quadrupole anisotropy 
causes CMB photons polarized by Thomson scattering with electrons at the last scattering surface 
{z ~ 1,100) and during the reionization epoch {z < 7) (Hu & White 1997). The amplitude 
of polarization is predicted to be 1 - 10% of that of the temperature anisotropies, depending on 
angular scales. The CMB polarization can provide useful information that is not much contained in 
the temperature anisotropy, such as the epoch of reionization or the tensor perturbations. Recently 
Keating et al. (2001) have reported the result of the POLAR experiment, giving upper limit of 
order 10 on the large-scale CMB polarization. The PIQUE experiment (Hedman et al. 2001, 
2002) have also obtained a similar upper limit at subdegree scales. De Oliveira-Costa et al. (2002) 
have tried to measure the cross-correlation Cj^ between temperature and i5^-mode polarization 
power spectra by cross-correlating the PIQUE and Saskatoon data (Netterficld et al. 1997). Many 
single dish and interferometry experiments such as MAXIPOL, POLAR, Polatron, COMPOSAR, 
CMB RoPE, DASI, and CBI are on-going to detect these faint polarized signals of CMB origin. 
Very recently, the DASI has reported a detection of the CMB E'-mode polarization and the TE 
cross-correlation by differencing the CMB fluctuations in two fields in 271 days of observation 
(Leitch et al. 2002b; Kovac et al. 2002). 

Since the Cambridge Anisotropy Telescope has detected the anisotropy on subdegree scales 
(CAT; O'Sullivan ct al. 1995; Scott et al. 1996), the Degree Angular Scale Interferometer (DASI; 
Leitch et al. 2002a; Halverson et al. 2002), the Cosmic Background Imager (CBI*^; Padin et al. 
2001; Mason et al. 2002; Pearson et al. 2002), and the Very Small Array (VSA; Taylor et al. 
2002; Scott et al. 2002) have also measured the CMB angular power spectrum down to subdegree 
scales. A desirable feature of the interferometer for CMB observation is that it directly measures 
the power spectrum, and that the polarimetry is a routine. In addition, many systematic problems 
that are inherent in single-dish experiments, such as the ground and near field atmospheric pickup 
and spurious polarization signals, can be significantly reduced. 

The observation of CMB polarization using close-packed interferometers arc also on-going or 
planned by the experiments like the DASI, CBI, and the forthcoming interim Array for Microwave 
Background Anisotropy (AMiBA^; Lo et al. 2001). They have full capabilities to probe the CMB 
temperature and polarization simultaneously. The feed horns are able to detect T, Q, U, and V 
Stokes parameters using complex linear or circular polarizers, aiming to detect CMB linear (Q and 
U) polarization at wave numbers 100 < ^ ^ 4, 000. 

This paper proposes the data analysis techniques for various strategies of the interferometric 
observation, especially when the ■ut'-space beam size is larger than the structures of the CMB power 
spectrum in the ^-space. An attempt is made to increase the resolution of the estimated power 
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spectrum in the ^-space. We have adopted three observational strategies, namely single pointing, 
mosaicking, and drift-scanning methods, and tested a few analysis methods extracting the angular 
power spectra from mock data for efficiency. 

The outline of this paper is as follows. In §2 and §3, beginning with a summary of CMB 
interferomctric observation, we describe a theoretical formalism for analyzing CMB intcrfcromctric 
data. A prescription for making mock CMB observations with interferomctric array is given. 
Application of each observational strategy is made in §4. A fast unbiased power spectrum estimator 
is introduced in §5. Finally, the summary and discussion are given in §6. 



Here we briefly summarize the CMB interferometry experiments (Hobson, Lasenby, & Jones 
1995; Hobson k Maisinger 2002, hereafter HM02; White et al. 1999a, hereafter W99). Only 
the observed CMB signal and instrument noise are considered, and other contaminations like the 
ground pickup or Galactic foregrounds are not taken into account. The sky is assumed to be flat 
since the field size surveyed by the interferomctric array is usually smaller than 10°. 

The visibility function of CMB temperature fluctuations sampled at a pointing position on 
the sky is the Fourier transform of the temperature anisotropy Ar(x) on the sky weighted by the 
primary beam ^(x — y^, i^), i.e.. 



(Hobson et al. 1995; W99). The factor dB^/dT = {2k%T'^ /cPh^) x^e^/(e^ - 1)^ is the conversion 
factor from temperature to intensity, where Bi,{T) is the Planck function, Ub is the Boltzmann's 
constant, h is the Planck constant, c is the speed of light, and x = hv/ksT. We assume T ~ 
TcMB = 2.725 lb 0.002 K (Mather et al. 1999). The baseline vector u is a variable conjugate to the 
coordinate x, and has a dimension of the inverse angle. In practice a vector connecting two centers 
of a pair of dishes determines a baseline vector u in unit of observing wavelength. For a small field, 
u = ^cot(7r/£) i/l-K (Hobson &; Magueijo 1996). Each frequency channel gives an independent 
set of visibility samples with different baseline vectors. 

The size of the primary beam A(x) determines the area of the sky that is viewed and hence 
the size of the map while the maximum spacing of the array determines the resolution, ^(x) is 
normalized to unity at the pointing center, A(0) = 1. The visibility in equation (1) can be also 
defined in the ui;-domain. 



2. CMB Interferometry Experiment 



2.1. Visibility 




(1) 
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where A{u,i') and AT(u) are Fourier transform pairs of A{-x.,v) and Ar(x), respectively. Let us 
define the generalized txv-space beam pattern Ay^(u, z/) at a pointing position from the phase 
reference point, 

Iy„(u,Z.)EEl(u,Z.)e*2.u.y^. (3) 

We adopt a beam pattern ^(u, u) given by W99, 



where Dx is the diameter of a dish (D) in unit of wavelength. It is defined by the autocorrelation 
of the electric field pattern of the flat-illuminatcd feed horn. This is equal to the autocorrelation of 
a pillbox of radius Dx/2 (van Waerbeke 2000). The primary beam pattern A(x.) obtained by the 
inverse Fourier transform of equation (4) is well approximated by a Gaussian function (for instance, 
FWHM c± 20' for a dish with D = 60 cm and the center frequency = 95 GHz) . 

Visibilities of CMB Q and U Stokes parameters, V-^^ and Vy^, can be similarly defined by 
replacing AT(x) with <5(x) and i7(x), respectively. The Q and U Stokes parameters depend on 
the observer's coordinate system. Under a rotation of coordinate system by tp on the sky, Q ziz iU 
transforms as a spin-2 tensor, {Q±iU) e'^'^''^{Q±iU) (Zaldarriaga k Seljak 1997). It is useful to 
deal with rotationally invariant polarization quantities, known as E- and B-mode polarizations. We 
follow the small-scale approximation of Seljak (1997) and Zaldarriaga &: Seljak (1997) to describe 
CMB temperature and polarization fields. 



Ar(x) = j (fu AT (u)( 
g(x) = J (fu [Eiu) COS 2(/)u - S(u) sin 2(/)u]e~^2™"', 
U{x) = J (fu [E{u) sin 2^^ + B{vl) cos 2^^]e~''^'''^-'' , (5) 

where E{\i) and B{\i) are the Fourier transforms of the E- and 5-mode polarization fields, and 0u 
is the direction angle of u. 



2.2. Visibility Covariance Matrix 

We are interested in measuring the angular power spectra of the CMB temperature fluctuations 
{Cj'^ = i{£ + l)Cj^/27r), polarizations (Cf^ and Cf^), and their cross-correlation (Cf^). We can 
relate them with the power spectrum S (u) defined in the uv-plane under the flat-sky approximation 



by, for example, £{£ + l)Cj^ /2it 2ttu'^Ste{u) where Ste{u) is defined by (aT{u)E*{v) 




Ste{u)S{u — v), and likewise for other components. 
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To estimate CMB temperature and polarization power spectra from the visibility data using the 
maximum likelihood analysis, it is essential to construct visibility covariance matrices. Combining 
equations (2) and (5) gives the visibility covariance matrix, 

= Qj! j d^wAy^{ui - w, Vi)A*y^{uj - w, Uj)SxY{^), (6) 

where X and Y correspond to T, Q, or [/, and i (j) is a visibility data index at a pointing position 

Ym (yn)- All possible combinations of 5xy(w), defined by ^X(v)y*(w)^ = 5xy(w)(5(v — w), are 
listed in Tabic 1. We assume no correlation of temperature and £^-mode polarization with i?-mode, 
which is expected in common inflationary scenarios (Hu &; White 1997). Changing the sign of Uj 
in equation (6) gives another complex covariance matrix, 

= J d^wAy^ {ui - w, n)Ay^ {uj + w, iyj)SxY i-w). (7) 

In practice, it is more convenient to use real quantities rather than complex ones. Separating 
the complex visibilities into real (V^) and imaginary (V^) parts and combining and Nmn 

appropriately, we obtain the followings (HM02, Appendix A), 

{Vyl{vii,Ui)Vy^Juj,uj)) = (Re[M^„]+Re[iV^J)/2, 

{V^Jui,Ui)Vl{uj,uj)) = (Re[M^„]-Re[iV^J)/2, 

{Vylini,Ui)Vl{uj,uj)) = (-Im[M^„]+Im[iV^„])/2, 

{V^Jni,Ui)Vy'l{uj,uj)) = (Im[M^J+Im[Ar^„])/2. (8) 



2.3. Simulated Observations 

We use the CMBFAST power spectra Cf^, Cf^, and the cross-correlation Cf^ (Seljak & 
Zaldarriaga 1996; Zaldarriaga, Seljak, & Bertschinger 1998) of a flat ACDM cosmological model. 
The model parameters are Q, = 1, Qa = 0.7, h = 0.82, il^/j.^ = 0.03, and = 0.975, which was 
given by the joint MAXIMA-BOOMERANG data analysis (Jaffe et al. 2001). We assume that the 
5-mode vanishes (-B(u) = 0). To generate AT, Q and U fields on a patch of the sky we use the 
Fourier transform relations in equation (5). For a given £ we construct a 2 x 2 matrix M that has 
Cf^ and Cf^ as the diagonal and Cj^ as the off-diagonal elements. We then Cholesky-decompose 
this matrix, M = LL^, where L is a lower triangular matrix, and assign random fluctuations Ar(u) 
and E(u) in each uv-cell using the following relation, 



Ar(u) 

E{n) 




(9) 
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Here ri and are random complex numbers drawn from the Gaussian distribution with zero mean 
and unit variance (for definite normahzation of Eq. 9 in a finite field, sec Wu 1999). The results 
are inverse Fourier transformed (Eq. 5) to yield AT(x), Q{x), and U{x.). These fields are then 
multiplied by the primary beam and Fourier transformed to give a regular array of visibilities (Eq. 
1). An observation is simulated by sampling the regular array at the u points specified by the dish 
configuration and observation strategy. The bilinear interpolation is used in the sampling. 

The instrument noise is simulated by adding a random complex number to each visibility 
whose real and imaginary parts are drawn from a Gaussian distribution with the variance of the 
noise predicted in a real observation. We use the sensitivity per baseline per polarization defined 
as (Wrobel & Walker 1999; Ng 2001) 

1 2kBTsys 1 

= A ; , (10) 

VsVa Aphys V2Az/ 

where Ai^ is the bandwidth, Aphys is the physical area of an elemental aperture, Tgys is the system 
noise temperature, and rjs (rja) is the system (aperture) efficiency. We assume ijg = rja = 0.8. 
Equation (10) is the sensitivity of a simple interferometer with a single real output consisting of 
the product of the voltages from two antennae. For a complex correlator, the noise is statistically 
the same in each of the two channels and thus the real and imaginary correlator outputs, cr^ and 
ai respectively, are equal with ar = cri = CTb^ where ai, = si,/-y/f^ and Tqcc is the correlator 
accumulation time. If we use the simultaneous dual polarizer with noise level per polarization, 
the T, Q, and U visibilities will obey Gaussian statistics characterized by zero means and standard 
deviations, ax = (tq = ajj = ab/V2. 

The instrumental noise of a single dish experiment has a strong correlation with itself, for 
instance due to the l/f noise statistics. But the noise of an interferometric array is randomly 
distributed over the baselines and observations. This is because the noise in each correlator output 
of a pair of dishes has independent random noise statistics. Therefore, the noise covariance matrix 
of the visibility data can be assumed to be diagonal, with the non-zero element simply given by 
the variance of the noise of each baseline component. 



3. Data Analysis 

3.1. Power Spectrum Estimation 

We summarize the basic scheme for the quadratic estimator of the CMB power spectrum based 
on the maximum likelihood analysis (Bond, Jaffe, &; Knox 1998; Tegmark 1997b; Tegmark & de 
Oliveira-Costa 2001). 

We construct the data vector V = (V^, V*3, V'^), with each visibility data point decomposed 
into real and imaginary parts. The observed visibility is the sum of contributions from the CMB 
signal and the instrument noise, Vi = Vsi + Ui {i = 1,2, ... , Np), where Np is the number of data 
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points including all the real and imaginary components of T, and U visibilities. The covariance 
matrix of the visibility data becomes 

Cij = {ViVj) = {VsiVsj) + {runj) = S^j + Nij. (11) 

The Sij and Nij are called signal and noise covariance matrices, respectively, and the CMB signal 
is assumed to be uncorrelated with the instrument noise. The matrix element Sij is obtained from 

equation (8), where we perform numerical integrations in equations (6) and (7). Assuming that the 
CMB T, Q, and U fields are multivariate Gaussian random variables, we get a maximum likelihood 
function, 

^«^'» = (2^^^^"="' (-V-C-V/2) (12) 

for a given set of parameters {Cb} = {Cf^, . . . , C^^g^-^, . . . , . . . , C^f }• Here the Cf, = 

27ru'^S{ui)) is an estimate of band power for a band centered on (6 = 1,2,..., N^), and is the 
total number of measured band powers. The superscript T denotes transpose of a vector. By using 
the Newton-Raphson method, after several iterations we can find a set of band powers that most 
likely fit the simulated data of the assumed cosmological model and thus maximize the likelihood 
function. The full quadratic estimator is 



1 



2 ^(F-^)66'TV 
b' 



(13) 



where the Nf, x A^^ matrix F is the Fisher information matrix defined as 

The square root of the diagonal component of the inverse Fisher matrix, (F~^)^^^, is the minimum 
standard deviation of the measured band power C5. The derivative of the signal covariance matrix 
with respect to a band power estimate, dS/dCh, is obtained by combining dMmn/dCh and dNmn/dCf, 
as in equation (8). For a band power of temperature anisotropy (C^^) with a band width ranging 
from |ub^| to lufejl, for example, 

and likewise for dN^n/dC^^ . Cosine or sine factors are included in the integrand for band powers 
of polarization and cross-correlation {C^^ and Cf^). 

For comparing a set of measured band powers {Ch} with cosmological models, we require a route 
to get band power expectation values (C^) from the model power spectrum Q. The expectation 
value for a band power spectrum estimate Cb is given by 

n) = E(F"')(.6' E CeFi,,=Y,{w^/^)Ce, (16) 
6' ^€aii, e'ev e 
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where F is the band-power Fisher matrix, and the F^, is a Fisher matrix element with bands of 
the individual multipoles. W^/£ is called the band power window function (Knox 1999; Pryke et 
al. 2002), and defined as 

^ = E w E ' E ^'/^ = 1- (17) 

This function can be considered as a filter that averages out the power spectrum with variance 
weighting scheme (see Tegmark 1997b; Tegmark & de Oliveira-Costa 2001 for a related discussion). 

Although the quadratic estimator is much faster than the direct evaluation of the likelihood 
function in finding the maxima, it still demands intensive computing power. The most time con- 
suming part is the calculation of F, which is 0{Np) operation. The visibility data obtained from 
mosaicking or drift-scanning usually involves a large amount of visibilities, but many of them are 
strongly correlated. One way to increase the efficiency of the calculation is to work with a subset of 
the data that contains most of the signal. This transformation can be obtained using the optimal 
subspace filtering, also known as the signal-to-noise eigenmode transform (see Tegmark, Taylor, & 
Heavens 1997 for a review). In this paper we follow the transformation method given by Bond et 
al. (1998, Appendix A). We first perform a whitening transformation, 

N — > N-V2p^N-^/2 = I, 

S — > N-^/2gN-i/2^ ^^g) 
V — > N-'^/^Y, 

where the whitening filter N~^/^ is the inverse of the Cholesky decomposition of N. Since N is 
diagonal, the above quantities reduce to simple forms, Sij, Sij/aiaj, and Vi/ai, respectively. By 
diagonalizing the signal covariance matrix using the similarity transformation, we get 

S R^(N-V2gN-i/2)R = ^ = diag(efc), 

V — > R^(N-^/2v) = v^ (19) 

where diagonal components of e are eigenvalues of the whitened signal covariance matrix, and R 
is an unitary matrix of which the A;th column vector is an eigenvector corresponding to the kth 
eigenvalue e^. Transformation of dS/dCb into the eigenbasis can be also performed in a similar 
way. In this new basis, the data T^^'s are uncorrelated variables with variances {V^"^) = I + Sk- 
The £fe's are usually called the signal-to- noise eigenmodes. In our analysis, we can keep only the 
top 20 - 30% of the modes and treat them as a new data set. This is an excellent approximation 
to the original dataset because the contribution of eigenmodes with lower signal-to-noise ratios to 
the CMB signal is relatively negligible. By replacing the covariance matrix and the data vector 
in equations (13) and (14) with the transformed ones in equation (19), we get a new version of 
quadratic estimator in this eigenmode basis (see Eqs. AlO and All of Bond et al. 1998). Although 
in this formalism the transformation into the signal-to-noise basis, which is 0{N^) operation, is 
the most expensive part, the remaining procedures do not require large computation. 
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3.2. Data Compression: Pixelization in real- and ttu-spaces 

It is essential to compress the visibility data, since the quadratic estimator for power spectrum 
requires very large operation proportional to 0{Np). The optimal subspace filtering described 
in §3.1 is one method for data compression. For a completely close-packed array, we can reduce 
the length of the data set by averaging all the visibilities with the same baseline vector for each 
component (T, Q, U) and observing frequency. For example, the 7-element AMiBA in hexag- 
onal configuration gives 21 baselines for each frequency channel, and they can be reduced to 9 
independent baselines (Fig. la). 

In many CMB experiments, a telescope scans the sky continuously many times over a specified 
survey area, and generates the time-ordered data stream. Analyzing this raw data directly to 
estimate cosmological parameters is intractable. The same problem arises in the drift-scanning 
mode of the CMB interferometric observation (see §4.3). Furthermore, if the platform of the 
interferometer rotates with respect to the survey field during the observation we will get a even 
larger visibility data set from a denser coverage of the iit;-plane with the moving baseline vectors. 
In this case pixelization in the real- and/or uv-spaces can serve as a data compression method (Wu 
2002). IIM02 have developed a pixelization scheme in the -plane (see also Myers et al. 2002). 
This can be understood as a commonly-used map-making process (Tegmark 1997a; Stompor et 
al. 2002). Here we generalize the uv-spsjce pixelization method of HM02 to include the real-space 
pixelization as well. 

The time-ordered visibility data obtained by an interferometric array can be modeled by the 
pixelized visibility data vector V^'^ multiplied by a pointing matrix A plus a time-ordered instru- 
mental noise vector n, 

ytod ^ ^ypix ^ ^ (20) 

In component notation it is 

Kt) ("^) = E ^W{p0^yf + ("^)' (21) 

where y(tj) and Uj^ are the observing position and the baseline vector at time tj, and yp and 
are the pixelized real- and ttf-space positions, respectively. The pointing matrix A is determined 
entirely by the observational or scanning scheme. Suppose that the total number of time-ordered 
visibility data at a frequency v is iVtod, and the length of the pixelized visibility is ATpjx. Then the 
-^tod X -^pix pointing matrix A is defined with matrix elements = 1 if (y(tj), u^) lies in the 

cell specified by (yp,u/), otherwise = 0. 

The noise covariance matrix of the time-ordered visibility data can be taken to be diagonal 
(see §2.3), 

(Nt)(ifc)(i'fc') = {f^y(ti){'^k)ny(t.,){nk')) = (TikSii'Skk'- (22) 
A solution VP'^ of equation (20) can be obtained by applying appropriate filtering matrix W, i.e., 

ypi^ = WV*°'^. (23) 
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We adopt the well-known minimum- variance filter (Janssen &; Gulkis 1992; Tegmark 1997a), 

W = (A^N^^A)"^ A^Nt~\ (24) 

This filter has an attractive property that the reconstruction error, 5Ve = VP'''— VP''' = (A^N^^ A)~^ 
A-^N^^n, becomes independent of VP'''. The noise covariance matrix in the pixelized visibility data 
becomes 

Np^<5Ve<5Vj) = (A^N,-iA)-\ (25) 

Generally, the reconstructed visibility VP'" is the minimum- variance estimate of VP'", and 
becomes maximum-likelihood estimate provided that the probability distribution of n is Gaussian. 
Since the time-ordered noise covariance matrix Nt is diagonal, the matrix equations for the pixelized 
visibility and its noise covariance (Eqs. 23 and 25) are simple bin-averaging processes, i.e., 

and 



(Np)(pO(p'/') = {l/ E '^/<'fkYpp'^ii'- (27) 

From equations (23) and (24) it can be understood that when W acts on V*"*^, the first operation 
A'^N^-'^V*"'^ maps the raw data into the pixelized visibility space; this step averages all data with 
coordinates (y(<i), u/t) that belong to the cell (yp, u;), weighed with noise variance cr|^. The second 
one, (A-^N^''^A) ^, is simply the normalization factor. 

The remaining issue is to determine the size of the pixel (both in real- and uw-spaces) before 
processing equations (26) and (27). In virtue of the sampling theorem the size of the ■ui'-pixel, 
Aupix, must not be bigger than a half of the FWHM of the beam, i.e., Aupix < AMf„hm/2 (Hobson 
et al. 1995; HM02). For single pointing observation, the half- width of the beam defined in equation 
(4) is given by 

For example, A'Ufwhm/2 ^ 36 for the CBI {D = 90 cm, vq ~ 30 GHz) and Aufwhm/2 ^ 77 for the 
AMiBA {D = 60 cm, = 95 GHz; Fig. 16). We can set a similar limit for the real-space pixel 
size, i.e., < A6'fwhm/2. 



4. Observational Strategies and Applications: CBI and AMiBA Examples 

In this section we describe three observational strategies and give simulation examples for the 
CBI and the AMiBA. The forthcoming interim AMiBA is an interferometric array of 7 elements. 
The characteristics of the AMiBA is summarized in Figure 1. Previous theoretical works and 
simulations on the AMiBA project can be found in Ng (2001, 2002), Pen et al. (2002), and Park 
& Park (2002). 
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4.1. Single Pointing Observation 

In CMB interferometric observation, the observational strategy depends on the characteristics 
of the interferometric array, especiaUy on the dish configuration and on the primary beam size in 
the -uv-plane. In the single pointing strategy, the array observes many points on the sky. This 
strategy is appropriate for the interferometer that has uniform u?;-coverage and sufficiently small 
primary uv-heam size. This interferometer like DASI can cover a wide range of angular scales and 
resolve the structures in the CMB power spectrum. One can increase the signal-to-noise ratio of 
visibility and reduce the sample variance of the CMB fluctuations by increasing integration time 
per field {tf) and the number of independent fields (nf), respectively (see Halverson et al. 2002 for 
the recent results of the single pointing observation of DASI) . 

In Figure 2a we plot the CMB power spectra of temperature anisotropy measured from the 
CBI mock survey adopting the single pointing strategy (nf = 30, t{ = 72 hours) with band power 
widths of Au = Aitfwhm/2 = 36 (filled circles) and 72 (open circles). We assume that 13 elemental 
apertures {D = 90 cm) are located on a platform with DASI configuration given by White et al. 
(1999b, Table 1). We also assume that CBI complex correlators (Tgys = 30 K) have 5 channels 
around 26 - 36 GHz. During the observation of a specified field, the platform does not change 
its parallactic angle so that all baseline vectors keep constant orientations relative to the field. To 
estimate band powers from the mock visibility samples, we used the quadratic estimator described 
in §3.1. At high i > 2,000, where noise dominates, strong band power anti-correlations exist. 
We define the band power correlation by / \J M^Myy , where Mi,y is the inverse Fisher matrix 
element (F^-*^);,;,'. Correlations between neighboring bands at lower I < 2,000 are from —57% to 
—35% for An = 36, while they are between —19% and —16% for Au = 72. 

We can reduce the correlations among the band power estimates by sampling the visibilities 
in the uv-\Aa,ne in such a way that the nu-space is better covered both in radial and transverse 
directions and the window functions of various baselines are densely overlapping. As a result of this, 
the resolution of the power spectrum in the ^-space (hereafter A^-resolution) can be increased. To 
demonstrate this idea, we consider a case with the CBI array whose platform rotates by 60° with 
steps of 10°, 6°, and 2° during the observation of a single field to get visibility samples with 6, 10, 
and 30 different orientations, respectively. These situations are very similar to the CBI deep fields 
case (Mason et al. 2002). Furthermore, the visibility data set with 6 orientations approximately 
equals to the data set in each pointing in the recent CBI mosaic observation (Pearson et al. 2002) 
that has 78 baselines, 10 channels, and 42 mosaic pointings. Here the total number of visibilities 
is N„ ^ 200,000, and the effective number of orientations is 200,000/780/42 ~ 6. We pixelize 
the uf-plane with Aup[x = 12, and use equations (26) and (27) to reduce the data size. Figure 
26 shows the power spectra measured from the pixelizcd CBI mock visibilities with 6 (nf = 15, 
= 24^ X 6; open triangles), 10 (nf = 9, tf = 24'^ x 10; filled circles), and 30 (nf = 6, tf = 12^ x 30; 
open circles) orientations. The anti-correlations between the power estimates at £ < 2,000 are 
significantly reduced as the number of orientations increases. Table 2 summarizes the band power 
correlations and the band power uncertainties in the power estimations of the above four cases. 
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To further investigate the above method of increasing A^-resolution via the uv-mosaicking, we 
have simulated mock observations using a two-clcment interferometer with D = 20 cm (Aufwiun = 
51.2) and two frequency channels (20 GHz whole bandwidth centered at uq = 95 GHz). The 
dishes change their separation from 40 cm to 60 cm, during which the platform is rotated to obtain 
nearly uniform visibility samples. We adopt different noise levels depending on the amplitudes 
of the model power spectrum, i.e., ai, = 4:3{Ci/Ceoo)^^^ mJy per visibility. Figure 3 shows the 
distribution of visibilities in the uv-plane (upper right box) and the band powers estimated from 
30 mock observations. The band width of the intermediate five bands is An = 21 (Ai = 132), and 
the band power correlations are between —29% and —18%, with an average of —24%. 

Hobson &; Magueijo (1996) and Tegmark (1997c) derived useful formulae for the uncertainties 
of the measured power spectrum. Suppose that S{ub) denotes the power S{u) averaged over a 
band ranging from Ub-^ to centered on Uh- The rms error for this band power is given by 
V2 [S{ub) + S'^°^^^{ub)) , where the first part is due to the sample variance, and the second is due 
to the noise variance. If there are N{ub) independent eigenmodes in a band and if nf independent 
fields are observed, the error will drop by a factor of y^n{N{ub). 

Let us grid the uv-plane with cell size uq. In the single pointing strategy, the natural choice 
for the cell size is uo ~ l/x/HJ where is the solid angle of the primary beam. Then, the noise 
power per uv-ceU. is given by 



s 



a% = = (29) 

'^VIS 

where s is the sensitivity of the instrument in units of //K sec^/^ given by s = {dBy / dT)~^ Sb/^s-, 
nvis is the number of visibilities per uv-ce\\, and /jyis is the number of visibilities per unit area 
in the u^^-plane; n^ns = Pvis'i^o- Suppose that the sampled area in the half of the uv-plane is 
^0 = ^(^62 ~ ^61)/^' ^'^'^ i^ number of visibilities in this area {N^ = AoPvis)- Given 
the detector sensitivity s, uv-plane coverage Aq, and total observation time ttotj we can specify a 
quantity to keep constant, 

= ^ = 'P^, (30) 
where w can be considered as the noise weight per unit uu-area, and ttot = UiU. 

Prom equations (29) and (30) the noise power spectrum can be obtained, S^°^^^{ub) — o-%Uo = 
w~^n{^l/Nv. If Aeff{ub) is the effective area occupied by independent visibility components, 
then N{ub) ~ A^* = Acs{ub)0,s for sufficiently wide band. Therefore, the uncertainty limit for the 
band power estimation becomes 



^^w«V^sw('<"''"^i- '''' 



How many pointings do we need to make to minimize the band power uncertainty? With t^ot and 
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Ug fixed, we minimize AS{ui,) with respect to nf to get the best choice 

If nf is fixed instead, the optimal observation time for each pointing becomes t{ = s'^AQfll/NvS{ub). 
This corresponds to the Umit making the noise and the sample variance contributions equal. 

Let Tis consider an example using a two-element AMiBA with a single frequency channel 
{Ny = N^Ij = 1). We measure only one £'-mode band power with a wide band centered at uj, 2± 206, 
with c± 129 and u^^ — 283, given by the FWHM limit of the uw-beam pattern (Eq. 28). The 
fluctuation power of E'-polarization is expected to be C^^ ~ 20 fiK^ in our ACDM model. With 
^tot = 6 months, we obtain the optimal parameters, nf = 32 and tf = 135 hours, which give the band 
power uncertainty AC^ = 27:ulAS{uh) ^ 10 /jK^. Our mock observations with these parameters 
has confirmed the validity of equation (31). 

The time required for a 4cj detection of the CMB polarization using the AMiBA with 7 dishes 
is worth discussing. Unless we increase the number of independent modes N{uh) within the band, it 
would be impossible to get the 4cr limit by observing only one field even with infinite integration time 
due to the sample variance. Here we consider only the shortest baseline components {Ny = 12, 
= 3). From equation (31) with a condition AS{ui)) = S{uh)/4:, a 4cj detection of the E- 
polarization can be achieved in a minimum total observation time of ttot = 20 days when nf = 43. 



4.2. Mosaicking 

The resolution in the itv-space is limited by the area of the sky that is surveyed, which is equal 
to the size of the primary beam in a single pointing observation. By mosaicking several contiguous 
pointing observations, we can increase A£-resolution in the band power estimates. Mosaicking does 
not increase the u-range, but simply enhances the resolution by allowing us to follow more periods 
of a given wave (see W99 for a complete discussion). For close-packed interferometers such as the 
AMiBA and the CBI whose uv-heam size larger than the structure of the CMB power spectrum, 
mosaicking is essential to increase the A£-rcsolution. In particular, since the 7-elcmcnt AMiBA has 
baselines of only three different lengths (see Fig. Ic), the uv-coverage is so sparse that the band 
powers estimated from a single pointing observation will have large sample variances and strong 
anti-correlations (Ng 2001). 

The increase of A^-resolution by mosaicking can be demonstrated by using the effective uv- 
heam pattern 

-| -^mo -1 -^mo 

Ieff(u, ^) = E ^y^("' ^) = 1^ E ^("' (33) 

-'''mo ■ 1 -''mo ■ 1 

where N^o is the total number of contiguous pointings and yj is the zth pointing position on the 
sky. Figure 4 shows the profiles of the effective txv-beams resulted from 7- and 19-pointing mosaics 
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with separation 68^0 in a hexagonal configuration. As the separation SOuio increases, the effective 
width becomes narrower. However, for SO^ao > 20', complex correlations among the band power 
estimates are expected due to the sidelobe effects (e.g., sec the long dashed curve in Fig. 46). Figure 
5 shows an example of the CMB temperature and polarization power spectra estimates expected 
in the AMiBA 19-pointing hexagonal mosaic mock observation in a period of 6 months. Using the 
AMiBA configured as in the Figure 1, we have made 10 independent mosaics with SOmo = 15' and 
an integration time of 22.7 hours per pointing (nf = 10, t{ = 22^7 x 19). The quadratic estimator 
and the optimal subspace filtering are applied to estimate the band powers with Am = 90 (see 
§3.1). The band power window functions and the corresponding band power expectation values 
(Cf,) calculated from equations (16) and (17) are also shown. It shows that the measured band 
powers excellently match with the theoretical values. 

For the mosaicking strategy, the uncertainty limit of the band power estimation in equation 
(31) is modified to 



where Jlmo is the mosaicked area of each field. Here we have adopted a new uv-pixel size Uq ~ l/f^mo) 
and a new noise power per uv-cell af^ = s^f^sfimo/'^visif- The optimal mosaicked area can be ob- 
tained by minimizing AS{ub) with respect to r^mo; giving the condition, S{uh) = 'w~^nfQ,s^ino/Nv 
Given a total observation time itot = nitf, the best choice for the mosaicked area becomes 

^mo = ^^^S{ub) = l^^l. S{Ub). (35) 

As an example, let us consider the visibilities measured by the shortest baselines in the AMiBA 
polarization observation, as described in §4.1. With 7 dishes, the total number of visibility outputs 

is = 12 with 3 independent baselines {N^ = 3), and the number of independent modes within the 
wide band {ub ~ 206, u ~ 129 - 283) becomes N{ub) ~ NKQrao/^s) = ^cff(^i6)^mo- From equation 
(35), we obtain an optimal square-shaped mosaic with size of ^opt; assuming a total observation 
time ttot = 6 months, 

^opt = ^Uo « (36) 



Hence, from equations (34)-(36), the uncertainty of the band power estimate only depends on the 
total survey area nfilmo for a given total observation time ttot) and is expected to be ACf, ~ 2 
fiK^ for ttot = 6 months. Note that choosing a smaller band width will increase ACb', if there 
are Ub subbands within a certain wide band, the size of error bar for each subband in the wide 
band would be roughly ACbul^^. Increasing r^mo (or decreasing nf) will give better results since we 
can get more information from larger scales and hence increase the resolution of the band power 
width Ai > vr/^opt- Table 3 summarizes the optimal parameter choices for the single pointing and 
mosaicking strategies, given a total observation time itot = 6 months. 
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Figures 60, 66, and 6c show the CMB power spectra expected in an AMiBA experiment when 
it performs 12 x 12 mosaic {S9„io = 15') observations of 3° x 3° area over five fields (rif = 5) with 
6 hours of integration time per pointing. We have assumed two frequency channels and used the 
quadratic estimator as well as the signal-to-noise eigenmode analysis in our computation. Except 
for the first and the last bands, the band power width is chosen to be Au = 45 {A£ = 283). The £- 
location of each band power is roughly estimated from the noise-weighted window function defined 
as 



where the prewhitening transformation in §3.1 is applied, and is the number of data points that 
contribute to the sensitivity at £ (Bond et al. 1998). Note that equation (37) is an approximate 
window function, and the Fisher matrix-derived window function (Eq. 17) should be used in 
general. The band power correlations are sufficiently small with the average of about —10%. We 
have made two sets of measurement of the band powers with band centers shifted by the half of 
the band width with respect to each other. The ACDM model power spectra convolved with a 
box of band power width weighted by the window function in equation (37) are plotted to guide 
comparison between the measured powers and the theory (thin curves and squares). Since there 
are about 4 bands {uh ~ 4) within the range of a wide band in this case {u ~ 129 - 283 centered 
at Ub ~ 206), one would expect the estimated E'-polarization band powers near Ub = 206 to be 
(7(Ch) ~ ACbnU"^ ~ 4 /iK^, which coincides with the uncertainty of the fourth band power in Figure 
66 (4.1 /xK^) and validates equation (34). 



In the mosaic strategy, one has to track every pointing of a specified mosaic field for a suf- 
ficiently long time. This inevitably suffers from severe ground pickups since the ground emission 
usually changes abruptly as the platform of the array turns around during its tracking. Drift- 
scanning is the simplest way to remove the ground contamination because an interferometer is 
insensitive to the ground emission that is a DC signal (Pen et al. 2002). In this strategy the direc- 
tion of the array is fixed while the sky drifts along the constant declination as the earth rotates. 
After some time, the array slews to the original starting position, and observes the sky along the 
same scanning path. In the drift-scan mode, the visibility is a function of time along a scan path. 

To estimate the power spectrum by using reasonable computer resources, the time-ordered 
visibility samples must be compressed (see §3.2). During the compression, the data experiences 
effects of smoothing over the pixel scale. If we consider only the real-space pixelization, a pixelized 
visibility T^(ufe) is related with the time-ordered visibilities T^^^(ujfc) by 




(37) 



4.3. Drift-Scanning 
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dT 
dT 



J (fvA{w)Wp{^^)e'^''^■y^AT{uk - v), (38) 



where Wp(v) is the Fourier transform of the smoothing filter Wp{y), y = y{t) is the drift-scan 
path, and denotes the real-space pixel position on the sky. 

We have made an idealized mock AMiBA experiment in which the array drift-scans a small 
square field of 2° x 2° size near the celestial equator (for flat sky approximation). The platform 
does not change its parallactic angle during the observation. This enables us to consider only the 
real-space pixelization in the data compression. After scanning a 2° path, the next path 10' apart 
from the previous one is scanned. The T, Q, and U visibilities are functions of time determined by 
the rotating speed of the sky, and are accumulated every 8 seconds. Also, according to equation 
(36), the AMiBA drift-scans a square field of 2° x 2° many times in a way described above for 
t{ = 15 days, and finally observes totally nf = 12 independent fields in 6 months. The time-ordered 
visibility samples are then compressed as described in §3.2 to the pixelized visibility data set with 
12 X 12 mosaic format and with 10' pixel size. The pixel smoothing effect is significant only in the 
scan direction. In this situation l^p(v) in equation (38) is top-hat smoothing filter 

where A6p = 10' is the real-space pixel size used in the data compression, and x denotes the scan 
direction. The visibility covariance matrix for the top-hat filtered drift-scan samples is obtained by 
modifying equations (6) and (7), 

^'rln = -Q^-Q^ J d'wAy^{Ui-W,l^i)A;^{uj-W,l^j)Wp{Ui^-W^)Wp{Uj^-W^)S^^ (40) 

and likewise for Nmn- The band power estimates from the visibilities obtained by the mock drift- 
scanning are shown in Figures 6d, 6e and 6/, with a band width Au = 62.5 except for the first and 
the last bands. The band power correlations between neighboring bands are from —12% to —10%. 



5. Fast Unbiased Estimator 

The AMiBA is expected to observe the sky over 100 deg^ with the drift-scanning strategy. Even 
if the pixelization is applied to compress the time-ordered visibility samples, it is still a formidable 
task for a data pipeline using the maximum likelihood quadratic estimator to compute the power 
spectra. In this section we propose a fast power spectrum estimator that can be implemented 
without recourse to large computing resource. The basic idea comes from the correlation function 
analysis, first developed by Szapudi et al. (2001), and from the MASTER method of Hivon et 
al. (2002). The visibilities are, by definition, the beam-convolved CMB Fourier modes (e.g.. 
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F^(u) oc J (fvA{v)AT{u — -v)). Therefore, the correlation between a pair of visibihties are directly 
related with the power spectrum (Eqs. 6 and 7). We introduce an approximate expression, 

ViVj - Nij ^ Sij = Bi,{b)Cb, (41) 

b 

where Bij{b) = dSij/dCf, {i, j = 1,2, - ■ ■ , Np-, b = 1, 2, • • • , TV;,). The ensemble average of the above 
equation restores equation {11), (ViVj) = Sij + Nij. 

First, we perform the whitening transformation of the noise covariance matrix N, the instru- 
ment filter function Bij{h), and the visibility data vector V as shown in equation (18). Since N is 
assumed to be diagonal, the transformed quantities are 5ij, Bij/aiGj (= -B^), and Vi/ai (= V^), 
respectively. Equation (41) then becomes 

Vrv;" - = Y,BZ{b)C,. (42) 

b 

For a given band, complex pair weighting might be required for each z_7-pair because of the uneven 
ui;-coverage. However, we adopt an uniform weighting scheme with equal weights. 

After a careful rearrangement of all components of V^V^" — 5ij, Bfj{b), and Cb into simpler 
forms, {vk}, {f3kb}, and {c^}, respectively, by merging an ij-pair to an index k, we can rewrite 
equation (42) as v = (3c. The fast unbiased estimator (FUE) of the CMB power spectrum is 
obtained as 

c= [/3^(3]~'l3^v, (43) 

where 

(/3^/3),,, = ^i?-(6)i?-(6'), (44) 

and 

iP^v)^ = ^ Br,{b) {Vrv^^ - d,,) . (45) 

The computation of equations (44) and (45), which are 0{Np) operations, is straightforward with- 
out need of any large memory and disk space. The FUE method is fast especially when the noise 
covariance matrix N is diagonal. Even for the non-diagonal noise covariance matrix, however, the 
FUE is still fast since the prewhitcning transformation is required just once. If the visibility samples 
are sufficiently separated so that they are considered as almost independent, the signal covariance 
matrix S (also Bij{b)) can be approximated to be diagonal. This approach, which is very simi- 
lar to a multi-band generalization of the Boughn-Cottingham statistic (Boughn et al. 1992), can 
accelerate the FUE speed from 0{Np) to 0{Np) operations. 

To constrain cosmological models using the FUE estimated power spectrum, we need to know 
the uncertainty for each Cb and the covariance between the band powers or the band power window 
functions. The uncertainty of each band power estimated by FUE can be obtained by a Monte 
Carlo simulation (see, e.g., Hivon et al. 2002). By fitting or interpolating the band power estimates 
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{Cf,}, we can obtain a smooth CMB power spectrum Cf^, from which many mock observations are 

made by inchiding CMB signals, instrumental noises, and other characteristics of the experiment. 
The mock data sets arc analyzed in the same way as the real data, giving rise to a set of power 
spectrum estimates {{C™^}, {C™^}, . . .}. The covariance matrix of the band powers 

C,^ ^ ( (Cf - (C,-)„„,,) (Cf; - (C,7)^„,,) )^^^^ , (46) 

is calculated. The uncertainty for each band power is then given by the square root of the 

1/2 

diagonal components of C, i.e., ACb = C^^ . 

We apply the FUE method to the AMiBA 19-pointing mosaic data to measure the band powers 
{Ch} with error bars estimated from 30 simulation data sets (open stars in Fig. 5). The band power 
estimates by FUE are quite consistent with those obtained by the quadratic estimator. However, 
most of the FUE error bars are bigger than those of the quadratic estimator. Note that the FUE 
is a sub-optimal quadratic estimator with compared to (S + N)"^ of the optimal quadratic 
estimator. 



6. Discussion 

We have simulated interferometric observation of CMB temperature and polarization fluctua- 
tions. For each observational strategy the data pipelines from the time-ordered raw visibility sam- 
ples to the CMB angular power spectra (C^^, C^^, and C^^) have been developed. The pipelines 
are composed of making mock observation, data compression, and power spectrum estimation. 
Data compression is achieved by pixelization of time-ordered visibilities in real- and uv-spaces by 
means of a common map-making process. This method can be applied to any kind of interfero- 
metric observation (sec §3.2). In estimating the band powers from the mock visibility samples, the 
optimal subspace filtering or signal-to-noise eigenmode analysis along with the quadratic estimator 
was used. By discarding the modes with low signal-to-noise ratios, we were able to reduce the data 
set to a manageable size. One drawback of the optimal subspace filtering is that while it conserves 
the information with signal-to-noise ratio higher than the limit of eigenvalue threshold, some useful 
information may be lost in certain cases. For instance, in measuring the CMB polarization power 
spectrum, if the band width is too small to keep sufficient amount of the signal compared with the 
noise level, the weak signal can disappear during the optimal subspace filtering. Therefore, we need 
to choose a wider band width to obtain a higher signal-to-noise ratio, especially at higher i region. 

The measured band powers are found to be quite consistent with the band power expectation 
values (Cf,) for the AMiBA 19-pointing mosaic (Fig. 5). This implies that our data pipelines 
are working reliably. Using the fact that the visibility contains direct information of CMB power 
spectrum, we have developed a fast unbiased estimator of the CMB power spectra (FUE, §5) that 
requires only 0{Np) operations. This method is very similar to the power spectrum estimation 
method using Gabor transform (Hansen, Gorski, &: Hivon 2002). The FUE also gives band power 
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estimates that are consistent with those from the quadratic estimator (see Fig. 5) . The FUE method 
does not require large computer resources. Given the precomputed quantities Bij{b) = dSij/dCf,, 
the computational speed is extremely fast. Even if the noise covariance matrix is highly non- 
diagonal, which is the usual case in real data analyses (e.g., handling constraint matrices to subtract 
the point source effect), the FUE method is still fast because the prewhitening transformation of 
Bij{b) and Vi is needed only once. 

Our main goal was to propose data analysis techniques for each observational strategy of a 
CMB interferometer, especially when the uv-beam size is larger than the scale of structures in the 
CMB power spectrum. Using the mock CBI single pointing observations, we have investigated the 
effect of rotation of the array platform on the band power correlations and the uncertainties of the 
band powers. Based on the results, summarized in Figure 2 and Table 2, we conclude that the 
band power anti-correlations can be reduced by rotating the platform and thus densely sampling 
the visibility plane. However, the uncertainties of the band power estimates slightly increase (when 
the total integration time is fixed). This is because the CMB signal is shared by the neighboring 
visibilities due to the finite beam size. In this way, wc can increase the resolution of the power 
spectrum in the ^-space down to a resolution limit A£ n/O given by the sampling theorem. 

Using the recent CBI result of single pointing observation, Mason et al. (2002) have shown a 
power spectrum with band width of A£ = 4-\/21n2/6'fwhm ~ 300. This limit for A£ is the FWHM 
of the visibility window function, which is proportional to the square of the Fourier transform of 
the primary beam with ^f^hm (Pearson et al. 2002). On the other hand, our choice for the band 
width is Ai = 226 {Au = A«fwhm/2 = 36). This is the limit given by the sampling theorem (A^ 
Tr/9). It is A£ = 41n2/^fwhm for a Gaussian primary beam, which is a factor of \/2 narrower than 
that adopted by CBI team. We show in Table 2 that a mock CBI observation with 30 different 
orientations results in about 20% anti-correlations between neighboring bands at £ < 1,000, and 
higher values at higher i > 1,000, while they arc 10 - 15% in Mason et al. (2002) due to the 
wider band width. The band power correlation at high £ regions can be reduced by more densely 
sampling visibilities with sufficient integration time. As shown in the example of u-u-mosaicking 
using a two-element interferometer (Fig. 3), the band widths of power spectrum can be reduced 
while keeping the band correlations at a tolerable level by increasing the number of rotation steps 
with increasing dish separation, and by assigning longer integration time to the visibilities at low 
CMB signal regions. For intermediate five bands, the band width and the average band power 
correlation are Ai = 132 (Au = 21) and —24%, respectively. The width is smaller than our 
resolution limit (A^ = 161, Au = Aufwhm/2 where Auf^hm = 51.2), also a factor of 1.7 narrower 
than the limit obtained by Pearson et al.'s formula {Ai ~ 230). 

The recent DASI power spectrum is measured from the single pointing observation without 
platform rotation (Halverson et al. 2002). The DASI band powers have the resolution of Ai 80 
(with 18 - 28% anti-correlations), which is broader than the resolution limit Ai = 4v^ln2/^fwhm = 

66 where ^fwhm = 3?4. We expect that the DASI single pointing observation with dense rotation of 
platform will allow higher resolution of about Ai = 50 at the similar level of anti-correlations. Since 
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the mosaicking is the most efficient method for increasing the resolution of the power spectrum, 

the combination of mosaicking and dense rotation of the platform followed by the uu-pixelization is 
thought to be the most ideal observational strategy for DASI- and CBI-type CMB interferometers. 

For each observational strategy, optimal parameter choices for the AMiBA experiment are 
discussed in §4, and summarized in Table 3. The 7-element AMiBA is expected to detect CMB 
polarization power spectrum near £ « 1,300 at 4cr level within 20 days by observing 43 fields. In 
AMiBA mosaicking with itot = 6 months, the optimal parameter sets are (^opt = 7°, nf = 1) or 
(^opt = 3°, nf = 5) for a minimum uncertainty of the ^^-polarization power spectrum. In fact, 
the optimal parameters strongly depend on the characteristics of the interferometer (e.g., 77^, rja, 
and Tgys) and on the ^^-polarization power spectrum C^^. Since we are considering the shortest 
baselines in deriving the parameters, the optimization is only for the sensitivity range of the shortest 
baselines {i < 2,000). At higher ^-range {i > 2,000) where the 7-element AMiBA has only a few 
baselines, the CMB polarized signal is expected to be very low. Therefore, we have chosen a wider 
band width for the last band for the 12 x 12 mosaicking and drift-scanning observations (see Fig. 6). 
To obtain a meaningful polarization power spectrum at high i region with narrow band widths, we 
need to increase the integration time or the number of baselines. This can be seen in the simulation 
of the 19-pointing mosaicking by AMiBA where the integration time per pointing is almost one day 
(see Fig. 5). Although the band widths are quite wide in the temperature power spectra in Figure 
6, we can measure temperature band powers independently with narrower band width because 
the signal-to-noise ratios of the T visibilities are very high, compared with those of polarization 
{A£ = 196, see open circles in Fig. 6c?). 

Among the three observational strategies that we have studied, the single pointing is useful 
for a detection of the CMB polarized signal while the mosaicking or the drift-scanning of a large 
area of the sky is essential for measuring the polarization power spectrum with high A£-resolution. 
The drift-scanning strategy is efficient for removing the ground contamination. It can also save 
half of the integration time when compared to the method of differencing two fields in the removal 
of the ground spillover adopted by the CBI experiment (Padin et al. 2001). In the drift-scanning, 
the survey region can have a shape for which the flat-sky approximation is inapplicable. Since the 
survey area drift-scanned by the AMiBA interferometer will be over 100 deg^, it is necessary to 
take into account the curvature of the sky. Our future work will deal with important issues such 
as the removal of Galactic foreground emission, the identification of radio point sources, and the 
subtraction of unresolved point sources. It is also important to study the topology of the CMB 
temperature and polarization fields to test the primordial fluctuations for Gaussianity (see, e.g.. 
Park & Park 2002). 
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Table 1. Ensemble Averages of the Conjugate Products of T, Q, and U Pairs in the u?;-Plane. 



XY <Jx(v)y*(w)^ = 5xy(w)(5(v - w) 



TT 


5tt(w) = Stt{w) 




TQ 


5r(5(w) = SxEiuj) cos 2</)w 




TU 


5tc/(w) = SxEiw) sin 2(/)w 




QQ 


'5qq(w) = See{w) cos^ 2(/)w + Sbb{w) 


sin^ 2^w 


QU 


Squ{'^) = See{w) cos 2^w sai2<j)^ - Sbb{w) 


sin 20W cos 20W 


UU 


Suu{^) = See{w) sin^ 2^w + Sbb{w) 


cos^ 2^w 
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Table 2. Band Power Correlations M^hi / y/MhbMh'b' and Relative Uncertainties ACb/Cb in the 

Mock CBI Observations? 



nf 




(6,6') 


(1,2) 


(2,3) 


(3,4) 


(4,5) 


(5,6) 


(6,7) 




Average 


30 


72^ 




-0.35 


-0.34 


-0.34 


-0.43 


-0.48 


-0.57 




-0.41 


15 


24h X gori. 




-0.25 


-0.29 


-0.27 


-0.40 


-0.44 


-0.43 




-0.35 


9 


24h X ;^o°"- 




-0.21 


-0.23 


-0.23 


-0.33 


-0.40 


-0.44 




-0.30 


6 


12'' X 30°"- 




-0.21 


-0.21 


-0.22 


-0.31 


-0.33 


-0.36 




-0.27 














ACb/Cb 










nf 


tf 


6 


1 


2 


3 


4 


5 


6 


7 


Average 


30 






0.15 


0.23 


0.11 


0.17 


0.14 


0.18 


0.21 


0.17 


15 


24h X gori. 




0.15 


0.27 


0.14 


0.17 


0.20 


0.18 


0.19 


0.19 


9 


24h X iQori. 




0.14 


0.28 


0.15 


0.20 


0.21 


0.22 


0.25 


0.21 


6 


12*^ X 30°"- 




0.13 


0.33 


0.20 


0.23 


0.24 


0.22 


0.25 


0.23 



^The total observation time is set to ttot = 90 days. 
^Mbb' is the inverse Fisher matrix element (F'"^)^^'. 

Note. — Only the first 7 band power estimates (An = 36 or A£ = 226) are considered. 
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Table 3. Optimal Parameters for AMiBA Single Pointing and Mosaicking Strategies? 



Single Pointing 




nf tf 




2 


32 135 hrs 


10 /iK^ 


7 


380 11.4 hrs 


2 /xK2 




Square Mosaicking (n^ 


= 7) 


^opt 


nf tf AUo 


iS9^o = 15') 


7° 


1 180 days 


28 X 28 


3° 


5 36 days 


12 X 12 



^The total observation time is set to ttot = 6 
months. 

^Number of dishes. 

^£^-mode band power uncertainty. The band 
width of the AMiBA shortest baselines is 
Au = 154, centered at uj, = 206. 
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Fig. 1. — (a) Close-packed hexagonal configuration of seven receivers (D = 60 cm) of the AMiBA 
on a single platform. Each feed horn can detect T, Q, U, and V Stokes parameters simultaneously 
using complex dual polarizer with Tgys = 70 K. Thin circles are boundaries of the physical dishes, 
and the shields between the 5 cm gaps of the adjacent dishes are drawn as thick circles. Due to 
the hexagonal configuration, the instrument gives 21 baselines with 3 different lengths in different 
directions, (b) The ■uw-coverage for a snapshot of the sky with the configuration shown in (a). Four 
frequency channels are assumed (t^o = 95 GHz, Au = 20 GHz whole bandwidth). The points are 
the ttt;-coordinates of the baseline vectors, and the circles have diameters equal to the FWHM of the 
uv-heam pattern, (c) The sensitivity of AMiBA at the center frequency vq = 95 GHz shown by the 
window functions for the three baselines with different lengths (Eq. 14 of W99). (d) The real-space 
primary beam pattern A(x) at the center frequency (solid curve), and a Gaussian function with 
FWHM ~ 20' (dotted curve) approximating the beam. 
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Fig. 2. — (a) The power spectra of the CMB temperature anisotropy measured from a mock CBI 
observation (nf = 30, = 72^) with band power widths Au = 36 (filled circles) and 72 (open 
circles), (b) Power spectra measured from mock CBI observations with 6 (wf = 15, t{ = 24^ x 6; 
open triangles), 10 (nf = 9, tf = 24^ x 10; filled circles), and 30 (nf = 6, tf = 12*^ x 30; open circles) 
orientations. Each figure is divided into two panels at low and high i regions. To avoid confusion, 
the open triangles and circles are slightly shifted to the left and right with respect to the filled 
circles, and open circles at ^ > 2000 are omitted. 
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Fig. 3. — The CMB temperature anisotropy power spectrum estimated from 30 mock observations 
by a two-element interferometer (D = 20 cm, two channels with Av = 20 GHz whole bandwidth 
centered at 95 GHz). The dish separation changes from 40 cm to 60 cm, and the platform is rotated 
to obtain uniform visibility samples. The distribution of visibilities in the uv-plane is shown in the 
upper right box. In the upper right corner of the box, the primary beam is shown as a circle 
with diameter Auf^hm = 51.2. Except for the first (Au = 25) and the last (Au = 50) bands, the 
band width of the intermediate five bands is chosen to be Au = 21 (Ai = 132). The band power 
correlations between neighboring bands are also given at the bottom of the figure. 
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Fig. 4. — Effective beam patterns expected (a) in a single pointing observation (solid curve) and 7- 
and 19-point hexagonal mosaicking (dotted and dashed curves) with pointing separation SOmo = 10', 
and (b) in 19-point hexagonal mosaics with pointing separations of S6„io = 10', 20', and 40'. The 
pattern of the hexagonal mosaic is similar to that of the AMiBA dishes shown in Figure la. Since 
the hexagonal configuration makes the imaginary components of the individual beam patterns along 
the Ux and cancel with one another, the effective beam patterns are real functions. 
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Fig. 5. — Band power estimates for (a) the temperature, {h) ^^-polarization, and (c) TE cross- 
correlation power spectra measured from the AMiBA 19-pointing hexagonal mosaic mock obser- 
vation (nf = 10, if = 22^^ 7 x 19, and 5^mo = 15')- Open circles with error bars are the results 
of the quadratic estimator (QE) and the signal-to-noise eigenmode analysis (S/N). For each band, 
the band power expectation value [Cb) (denoted by an open diamond) and the corresponding band 
power window function /J* = W^/l (denoted by a solid or a dotted curve in the bottom region) are 
also shown. The i'-location of each band power (£cfr) is found from 1^,^ = ^^f-f^ / f\ (Bond et 
al. 1998). The band power estimates obtained by the fast unbiased estimator (FUE), right-shifted 
with respect to l^s and denoted by open stars, are consistent with those of the quadratic estimator. 
The FUE error bars are estimated from 30 simulations (see §5). 
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Fig. 6. — (a) Temperature, (b) ^'-polarization, and (c) TE cross-correlation power spectra measured 
from an AMiBA 12 x 12 mosaicking observation with ^opt = 3°, nf = 5, and 56mo = 15' for 6 months. 
Band powers are measured at two sets of band centers (filled and open circles) shifted by a half 
of the band width with respect to each other. Band power estimates from the data obtained by 
the AMiBA drift-scanning strategy (^opt = 2°, nf = 12, and ttot = 6 months) followed by the 
12 X 12 pixelization with A^p = 10' are shown in {d), (e) and (/). Open circles in the panel (d) are 
the temperature power spectrum measured with the band width (Au = 31.25) twice smaller than 
that of the filled circles. The location of each band power is estimated using the window function 
defined in equation (37). Thick curves are the ACDM model power spectra. The model power 
spectra smoothed by a box of band power width {Au = 45 for mosaicking, and Au = 62.5 for 
drift-scanning) weighted by the window function are also shown (thin curves, and filled and open 
squares) . 



